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Probing the geometry of the universe is one of the most important endevours in cosmology. 
Current observational data from the Cosmic Microwave Background anisotropy (CMB), galaxy 
surveys and type la supernovae (SNe la) strongly constrain the curvature of the universe to be close 
to zero for a universe dominated by a cosmological constant or dark energy with a constant equation 
of state. Here we investigate the role of cosmic priors on deriving these tight bounds on geometry, 
by considering a landscape motivated scenario with an oscillating curvature term. We perform 
a likelihood analysis of current data under such a model of non-trivial geometry and find that 
the uncertainties on curvature, and correspondingly on parameters of the matter and dark energy 
sectors, are larger. Future dark energy experiments together with CMB data from experiments 
like Planck could dramatically improve our ability to constrain cosmic curvature under such models 
enabling us to probe possible imprints of quantum gravity. 



I. INTRODUCTION 

Exquisite measurements of the Cosmic Microwave 
Background (CMB) spectra, Large Scale Structure 
(LSS) and of the expansion history of the universe, 
by various experiments including_WMAP, SDSS, 2dF 
survey, and SNla surveys P, [2,H, Q, [B[ have pinned 
down the spatial curvature of the universe to K < 0.01 
0, 0] in a universe with a cosmological constant . 

Energy content and geometry both contribute to 
the Hubble expansion rate H. The three are so closely 
intertwined that no independent measurement of H 
versus il m , Q\, or n can yet be performed. In inter- 
preting the data for any of these components, we have 
to be aware of the prior assumptions in our models re- 
garding the other components. That is to say that our 
interpretation of data is not model independent. 

What conclusion can be made about the geometry 
of the universe, on the basis of current data? The 
answer to this question depends on the cosmological 
model and on the dependence of the constraints on the 
spatial geometry of the universe on prior assumptions 
regarding other relevant energy components. Here we 
explore this issue by considering a model of nontrivial 
geometry, where the curvature does not take a con- 
stant value, but rather is a function of time. In our 
model the curvature is given by an oscillating function 
with a Hubble time period. Such models can be moti- 
vated by the dissipative dynamics of the wavefunction 
of the universe on its classical path on the background 
of the landscape of string theory. 



II. MOTIVATING THE CLASS OF 
OSCILLATING CURVATURE MODELS 

The discovery of the acceleration of the universe has 
become one of the central themes of current investiga- 
tion in physics. Unfortunately due to the degeneracy 
among cosmic parameters, determining its nature and 
equation of state from astrophysical data depends cru- 
cially on our assumptions for the matter content and 
curvature of the universe. Recent analysis of SDSS0] 
and WMAP Q has been reported to indicate that the 
geometry of our universe is extremely close to flat to 
within 1% in an LCDM universe 

Here we would like to investigate the robustness of 
these conclusions about the curvature of the universe 
by presenting a highly nontrivial model where the cur- 
vature term is not a constant but a function of time, 
oscillating every Hubble time. This model is inspired 
and motivated from the proposal for a dynamic selec- 
tion of the initial conditions for our universe from the 
landscape phase space [8|, [9( as summarized below. 

In [H, Q we included the backreaction of superhori- 
zon massive perturbations on the initial wavepacket 
for our universe. Solving a Master Equation we stud- 
ied role of the backreaction term on the decoherence 
of our initial patch from the other WKB branches on 
the landscape. The time evolved nonlocal entangle- 
ment of our patch with others outside the horizon at 
late times [13, EU were then investigated in [l2| with 
the conclusion that some of those traces imprinted on 
CMB and LSS are within observational reach. Within 
this formalism, we now allow for an initial curvature 
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term in the Friedman equation, k = 1 and consider the 
effect that the backreaction term has on ciosed geome- 
tries. Backreaction shifts the energy of the wavepacket 
therefore its classical trajectory on the landscape. The 
shifting of the classical path for our wavefunction can 
be seen by integrating out the Master Equation pre- 
sented for this case as was described in [8, 1(| 12|. 



H{a) + H„ 



(1) 



where —H(a) = a(a? + k) is the hamiltonian cor- 
responding to gravitational degrees of freedom (ob- 
tained from the Einstein- Hilbcrt action), with a the 
scale factor, and k is the initial curvature for that 
classical trajectory. H m is the matter hamiltonian 
corresponding for example to the inflaton energy and, 
H n is the backreaction energy corresponding to the 
superhorizon matter perturbations labelled by the 
wavenumber n. 

To get a rough idea of the shifting of the trajec- 
tory, let us assume that e = a 3 [<p 2 + m 2 (j) 2 ] is a con- 
stant of motion and thus integrate out Eqn[TJ When 
the backreaction term is not included, EqnQ] in the 
case of closed universes K = 1, gives a turning point 
when a = at a = a n where a n n — e. The first 
term in the backreaction energy H n ~ n 2 /a 2 has the 
same dependence on the scale factor a as the cur- 
vature. Including the backreaction term when in- 
tegrating out Eqn[Tj results in a lower energy since 
H(a) — > H(a) — ^i? n , thereby shifting the classical 
trajectory of the wavepacket of our universe. This re- 
sult of the shifting of the classical trajectory of our 
universe's wavepacket by interaction with a field, (in 
this case the entanglement of gravitational and matter 
degrees of freedom through the term ^if n ), is well 
known in particle physics where the energy of a quan- 
tum particle gets shifted by interaction with a classical 
field which results in a shifts of the particle's trajec- 
tory and momenta. The details of the calculations for 
the strength of this interaction in the case of quan- 
tum cosmology can be found for example in o r 
when applied to the landscape of string theory in 12J . 
The result is that everytime a closed universe goes 
through its turning point given by EqnJT] by putting 
a = 0, that is every Hubble time, then it will emerge 
through a shifted trajectory, due to the correction of 
H(a) by the backreaction term described here. 

The modification in energy corresponding to this shift 
can be absorbed into the curvature term since the time 
dependence of H n is similar to that of the k term. 
Therefore, to local observers bound to our visible uni- 
verse, the effect of the shifting of the classical trajec- 
tory for our universe in the phase space, appears as an 
induced oscillating curvature with a period of Hubble 
time. Of course no observers would survive the emer- 
gence through the turning point in the cycles of the 



trajectory. However the reduced oscillation of the cur- 
vature in the previous cycle may leave its imprints on 
astrophysical observables of the current cycle, which 
we aim to study here. Based on the integration of 
EqnfTJ we expect the curvature term motivated by 
this scenario, to be a function of the total energy con- 
tent of the universe H m and oscillate, in each cycle 
between the two turning points in the trajectory, with 
a Hubble period. 

While there are different ways of phenomenolog- 
ically implementing the oscillating curvature model 
described above, there are constraints based on the 
considerations of observables in the model. In order 
to make predictions that can be compared with data, 
we need to compute the coordinate distance r(z) and 
the age of the universe t(z) as functions of redshift. 

These functions, r{z) and t(z), are found by consid- 
ering the radial, null geodesic of the Robertson- Walker 
metric. But now the curvature constant n is replaced 
by a function — x(a) that oscillates with cosmic time: 



ds 2 



dt 2 - a 2 (t) 



dr 2 



1 + x(a)r 2 



0. 



(2) 



where a(t) is the cosmic scale factor, and the scaled 
coordinate distance r = H$r / c. Hence we have 



,2 l + x(a)r 2 
(r ) = 



(3) 



a 4 E 2 {a) ' 

where / denotes differentiation with respect to a, and 



E 2 {a) = 



_ (H{z)\ _ n r n 



Hn 



x(a) 



= —r H r + . 

a 4 ar ar 



(4) 



For closed universes, r(z) is not a monotonic function 
of a. Differentiating Eq.J3]) with respect to a gives 

_„ x'r 2 + r' [2rx - r' (fi m + 4a 3 f2 A + x'a 2 + 2ax)] 
r = 2a 4 E 2 (a)r' ' 

(5) 

In a closed universe, the coordinate distance r reaches 
its maximum value at r' = 0. Note that r" is only 
finite at r' = if and only if x' = at r' = 0. This can 
only be satisfied if x ^constant (the usual constant 
curvature case), or x = x(r). 

Based on the dissipative dynamics of the shifted 
cycles of the universe described above, we consider 
the following heuristic model that captures the desired 
features and satisfies the above constraint 



x(a) 



{l + zf 



= n k -Asin(Br), 



(6) 



where J\- = ^fc(O) = 1 — fi r — O m — J7a, ^r, ^m> f^A 
denote the present day density fractions of radiation, 
matter, and vacuum energy, and r denotes the coordi- 
nate distance from the observer at z — to redshift z. 
A and B are dimensionless constants. Note that the 
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FIG. 1: Oscillating curvature models with B = 0.5, 1, and 
2. 



FIG. 2: The coordinate distance r(z) for oscillating cur- 
vature models in Fig.l, with the same line types. 



conventional model with constant curvature is recov- 
ered for A = 0, x(a) = tt k {z)/(l+z) 2 = fl k = -k/H$, 
where n is the curvature constant. 

Except for the special case of constant curvature 
(A = 0), r(z) is found by numerically solving the sec- 
ond order differential equation in Eq. ([5]) , with the ini- 
tial condition that ata=l,r = 0,r=— 1. 

Fig.l shows models with B — 0.5, 1, and 2 respec- 
tively. Fig. 2 shows r(z) for the models in Fig.l, with 
the same line types. 

The age of the universe is given by 

t(a)=H -iP^ r (7) 
Jo aE(a) 

once rife(z), which now depends on z through r(z), 
has been found numerically. 



III. OBSERVATIONAL CONSTRAINTS ON 
OSCILLATING CURVATURE AND THE 
ENERGY CONTENT 

We use current observational data to constrain the 
oscillating curvature model given by Eq.©. Follow- 
ing the approach of we assume the HST prior of 
H = 72 ± 8 (km/s)Mpc- 1 0L use 182 SNe la (from 
the HST/GOODS program [TH] , the first year Super- 
nova Legacy Survey [16j . and nearby SN la surveys) 
[HI, CMB data @, and the SDSS measurement of the 
baryon acoustic oscillation scale [ItJ ■ We use the CMB 
data in the form of the CMB shift parameters R = 



y/il m H^ r(zcMB) and l a = nr (zqmb ) ItAzcmb ) de- 
rived from WMAP three year data by [l3| • 

We run a Markov Chain Monte Carlo (MCMC) 18] 
to obtain O(10 6 ) samples for each set of results pre- 
sented in this paper. The chains are subsequently ap- 
propriately thinned. 

Due to the degeneracies between A, B, and f2&, (A, 
B, n k ) are not well constrained when they are all al- 
lowed to vary. To illustrate the effect of oscillating 
curvature, let us study the class of models given by 
Eq.© for fixed representative values of B, while al- 
lowing A and fl k to vary, along with fl m , flbh 2 , and 
h (see [13|). The parameters estimated from data are 
(Cl m , Qbh 2 , h, fik, A). It should be noted that A 
affects the overall amplitude of the curvature term, 
while B plays the role of its oscillating frequency. The 
case B — 1 would correspond to fl k oscillating every 
Hubble time. 

Figs|3][6] show the joint confidence contours in the 
plane (O m , fife) and (fife, A) for B = 0.5, 1, 2, and 5 
respectively. The inner and outer contours correspond 
to 68% and 95% confidence levels respectively. 

As can be seen from the plots given in Figs|3][6l 
current data allow models in which the curvature of 
the universe oscillates with cosmic time. The allowed 
range of the current curvature density ratio Ofc is sig- 
nificantly increased compared to the case of constant 
curvature. 

The bounds derived from the WMAP three year 
data and galaxy survey data from the SDSS 0] give for 
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FIG. 3: Joint confidence contour plots for (Q m , £lk) and 
(fifc, A) for B = 0.5. The inner and outer contours corre- 
spond to 68% and 95% confidence levels respectively. 




FIG. 4: Joint confidence contour plots for (Q m , fife) and 
(fjfc, A) for B = 1. The inner and outer contours corre- 
spond to 68% and 95% confidence levels respectively. 



the case of constant curvature, H k = —0.005 ± 0.006, 
(2dF data 0] also give similar results) @. Compar- 
ing these bounds to the case of oscillating curvature 
models, we find that the constraints on the geome- 
try of the universe change significantly, now we have 
n k = 0.097 ± 0.210 for B = 0.5, n k = -0.037 ± 0.092 
for B = 1, Q k = 0.014 ± 0.035 for B = 2, and 
ftfc = 0.000 ± 0.017 for B = 5. The constraints on 
the Ofc and A become more stringent as B increases. 
This is as expected, since B is the curvature oscilla- 
tion frequency. For large B, the cumulative effect of 
the oscillating curvature decreases. It is very interest- 
ing that when the period of the curvature oscillation 
B^ 1 becomes larger than a Hubble time, the range 
of the allowed values for fl k and fi m , at 95% confi- 
dence level in agreement with data, shows a drastic 



increase. An oscillation in the curvature with a pe- 
riod larger than the age of the universe, a case which 
locally would appear as nearly a constant while being 
globally notrivial, the time dependence of which would 
otherwise not be captured by data, does in fact con- 
tain a significant deviation from the priors of a simple 
LCDM model with constant or zero curvature. This 
is one of our important results: a highly nontrivial ge- 
ometry on scales larger than the horizon can lead to 
a very different interpretation of data. 

Let us now investigate the implications of the oscil- 
lating curvature for the dark energy equation of state 
w. This is done by conducting a likelihood analy- 
sis using MCMC of the oscillating curvature model, 
Eq.©, assuming a constant dark energy equation of 
state wx(z) = w. The parameters estimated from 
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FIG. 5: Joint confidence contour plots for (fl m , fifc) and 
(SI*:, j4) for B = 2. The inner and outer contours corre- 
spond to 68% and 95% confidence levels respectively. 



FIG. 6: Joint confidence contour plots for (fl m , fife) and 
(fife, j4) for B = 5. The inner and outer contours corre- 
spond to 68% and 95% confidence levels respectively. 



data are (f2 m , il&/i 2 , h, Slfc, A, w). Fig [7| shows the 
joint confidence contours of (f2 m , Qk) and (fifc, A), 
and FigO shows the joint confidence contour of (fife, 
w), for B = 5. As expected, adding w as an addi- 
tional parameter to be estimated from data notably 
increase the uncertainties on estimated parameters, 
especially (f2fc, A, w). For example, for B = 5, 
f2fc = 0.074±0.068 when w is included as an estimated 
parameter, compared with flk = 0.000 ±0.017 for set- 
ting w = — 1 (a cosmological constant). For larger 
values of the period -B -1 , the uncertainties on A, 
w) increase significantly. Notice that the bounds on 
w, as shown in Fig. 8 for the case B = 5 can be as 
large as —1.04 < w < —0.71 at 95% confidence level. 
These bounds on w should be contrasted to the con- 
straints derived in (lfjj where a prior of £lk = was 



assumed. Clearly, current precision cosmology data is 
not sufficient in pinning down the equation of state 
for dark energy when the geometry of the universe is 
nontrivial. 



IV. DISCUSSION 

We have studied constraints on the parameters of a 
landscape motivated cosmological model in which the 
curvature of the universe oscillates with cosmic time 
(see Eq.©). Such a model is motivated from the pro- 
posal for a dynamic selection of the initial conditions 
for our universe from the landscape phase space. Thus 
an analysis of the kind performed here could lead to 
the implicit detection of quantum gravity effects. 
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FIG. 7: Joint confidence contour plots for (fi m , f2fe) and 
(fifc, >1) for B = 5. The inner and outer contours corre- 
spond to 68% and 95% confidence levels respectively. 



We have used CMB data in the form of the shift 
parameters R and l a extracted from WMAP three 
year data by [l3j], together with the SDSS measure- 
ment of baryon acoustic oscillation scale [l7|, and 
SN la data from HST and ground-based observations 



[15l llq . From the bounds derived on the parame- 
ters of this model we find that currently a simple flat 
model, which is a special case of the above model, re- 
mains a good bet; such a conclusion will be supported 
further by model selection arguments [25[. Allowing 
for nontrivial geometry leads to greater uncertainties 
in our knowledge of the present day curvature and 
matter density ratios ttk and f2 m , as can be seen in 
FigslHHHJ- An oscillating curvature term also signifi- 
cantly changes the bounds on the dark energy equa- 
tion of state w, as seen in FigJSJ 
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FIG. 8: Joint confidence contour plots for (Qk, w) for 
B = 5. The inner and outer contours correspond to 68% 
and 95% confidence levels respectively. 



It would be interesting to look for the imprints of 
such a model as data get better. Future dark energy 
exp eriments from both ground and space [2ll . 22[ 23 , 
Q, together with CMB data from Planck [2(|, should 
dramatically improve our ability to constrain cosmic 
curvature, and probe possible imprints of quantum 
gravity. 
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